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Abstract 

We write a 3D equation for three fermions by combining the three two-body potentials obtained in 3D 
reductions (based on a series expansion around a relative-energy fixing "approximation" of the free 
propagators) of the corresponding two-fermion Bethe-Salpeter equations to equivalent 3D equations, 
putting the third fermion on its positive-energy mass shell. In this way, the cluster-separated limits 
are exact, and the Lorentz invariance / cluster separability requirement is automatically satisfied, 
provided no supplementary approximation, like the Born approximation, is made. The use of positive 
free-energy projectors in the chosen reductions of the two-fermion Bethe-Salpeter equations prevents 
continuum dissolution in our 3D three-fermion equation. The potentials are hermitian below the 
inelastic threshold and depend only slowly on the total three-fermion energy. The one high-mass 
limits are approximately exact. 

This "hand- made" three-fermion 3D equation is also obtained by starting with an approximation 
of the three-fermion Bethe-Salpeter equation, in which the three-body kernel is neglected and the two- 
body kernels approached by positive-energy instantaneous expressions, with the spectator fermion on 
the mass shell. The neglected terms are then transformed into corrections to the 3D equation, in 
three steps implying each a series expansion. The result is of course complicated, but the lowest-order 
contributions of these correction terms to the energy spectrum remain manageable. 

We also present some other 3D reduction procedures and compare them to our's: use of Sazdjian's 
covariant approximation of the free propagator, 3D reductions performed by a series expansion around 
instantaneous approximations of the kernels instead of "approximations" of the propagators, Gross' 
spectator model. 

PACS ll.10.Qr Rclativistic wave equations. 

PACS 11. 10. St Bound and unstable states; Bethe-Salpeter equations. 
PACS 12.20.Ds Specific calculations and limits of quantum electrodynamics. 
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1 Introduction 

In the treatment of the three-body problem, there is a large gap between nonrelativistic quantum mechan- 
ics (Schrodinger equation) and relativistic quantum field theory (Bethe-Salpeter equation |l|, || ) . Starting 
from the Schrodinger equation one can of course replace the free part of the hamiltonian by its relativistic 
form. Starting from the Bethe-Salpeter equation, one can try to eliminate the two relative time variables 
to get finally a Schrodinger equation with a compact potential plus a lot of correction terms of various 
origins. Instead of the Schrodinger equation, one can use the Faddeev equations. These equations can 
be derived from the Schrodinger equation, or, in a more general form, from the Bethe-Salpeter equation. 
Faddeev equations for the transition operator give the various scattering and reaction matrix elements 
and the poles (to be computed by basically nonperturbative methods) of this operator in the total energy 
give the spectrum of the three-body bound states. 

Schrodinger's equation is "relatively easy" to solve, but this zero-order approximation does not reflect 
several important properties and symmetries of the studied physical system. These could be recovered, 
in principle, by incorporating higher-order contributions (in general an infinity of them). 

As an intermediate step between nonrelativistic quantum mechanics and quantum field theory we 
shall search for a 3D equation built with the sum of the relativistic free hamiltonians plus two-body 
and perhaps also three-body interaction potentials. Such equations are closely related to the systems 
of coupled Dirac or Klein-Gordon equations of constraint theory [3-7], which are still 4D equations but 
exhibit a simplified dependence in the relative times, which could be completely eliminated to get a single 
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3D equation. Our 3D equation will be a strategic link in a chain of approximations between a manageable 
equation (such as Schrodinger's equation) and a full 3D equivalent (at least for what concerns the bound 
states spectrum) of the original Bethe-Salpeter equation. We shall try to keep a maximum of properties 
and symmetries of the original Bethe-Salpeter equation in our 3D equation. For the properties which 
we shall be unable to satisfy exactly in our 3D equation, we expect that the inclusion of higher-order 
correction terms will progressively improve the situation. We shall thus try to make our 3D equation 
satisfy the following list of requirements: 

— Correct nonrclativistic limit. 

— Lorentz invariance and cluster separability. It is always possible to render an equation Lorentz 
invariant by working in the general rest frame (center of mass reference frame) and by building invariants 
with the total 4-momentum vector, although the result may be unelegant and artificial. The Lorentz 
invariance requirement becomes a tool when combined with the cluster separability requirement: when 
all mutual interactions are "switched off", we must get a set of three free Dirac equations. This total 
separability can easily be obtained by using as hamiltonian the sum of three free Dirac hamiltonians and 
interaction terms. The real difficulty appears when only the interactions with fermion 3 (for example) 
are switched off. If we want a full cluster separability, the resulting equation for the (12) cluster can not 
refer to the global center of mass frame anymore, as the momentum of fermion 3 enters in the definition 
of this frame. 

Lorentz invariance and cluster separability can be explicit or implicit (via rearrangements) . The best 
known example of implicit Lorentz invariance is a free Dirac equation solved with respect to the energy: 
it becomes explicitly covariant by multiplication with the /3 matrix. Other implicit Lorentz invariances 
are not that trivial. For example, the 3D reductions of a Bethe-Salpeter equation are implicitly covariant, 
provided the scries generated by this reduction is not truncated. 

In fact, the homogeneous Bethe-Salpeter equation from which we shall start does not obey cluster 
separability: this equation is valid for three-fermion bound states only, i.e. for a total energy below all 
the continua. Our 3D equation will thus be used in the computation of three-fermion bound states. A 
priori, it will not necessarily give also the correct scattering amplitudes or obey cluster separability, even 
with all higher-order terms included. We shall nevertheless try to get a 3D equation which is also valid 
in the continuum region. 

— Hermiticity and total energy independence of the interaction terms. In the two-body problem, 
nonhcrmitian interaction terms can be hermitian with respect to a modified scalar product, or made 
hermitian via a rearrangement of the equation. In the three-body problem these rearrangements could be 
more complicated. The hermiticity and the independence on the total energy are linked features, as one 
of these is often achieved at the expense of the other one. Energy depending interaction terms destroy 
some of the advantages of the use of an hermitian hamiltonian, such as the mutual orthogonality of the 
solutions, and leads to modify the usual perturbation calculation methods. The 3D potentials deduced 
from field theory are generally energy dependent, at least in the higher-order terms. We shall require 
hermiticity and energy-independence (or slow energy dependence) below the inelastic threshold in the 
lowest order terms at least. 

— Correct heavy mass limits. When the mass of one of the fermions becomes infinite, its presence must 
be translated in the equations by a potential (Coulombian in QED) acting on the other fermions. For the 
two-body problem the "one-body limit" in QED is indeed a Dirac-Coulomb equation or a rearrangement 
of it (for example when a projection operator is introduced in order to avoid continuum dissolution - see 
below) . In the Dirac-Coulomb equation, the Coulomb potential is already given by the limit of the Born 
term, as the higher-order crossed and ladder terms cancel mutually at the one-body limit. In contrast, 
the one-body limits of the rearrangements contain contributions from all terms [8-12]. 

— Perturbative approach. It must be possible to start with a manageable approximate equation and 
indefinitely improve the approximation of the measurable quantities (with respect to the uncalculated 
predictions of the here assumed exact Bethe-Salpeter equation) by adding higher-order contributions. 

— Solution of the continuum dissolution problem. In the relativistic equations for several relativistic 
particles, the physical bound states are degenerate with a continuum of states combining asymptotically 
free particles with opposite energy signs. This often neglected fact forbids the building of normalizable 
solutions in the N > 2-body problem (including the two-body plus potential problem. In the pure two- 
body problem the mixing is prevented by the conservation of the total momentum). The usual solution 
consists in including positive-energy projection operators into the zero-order propagator [13-16]. The 
modified equations must of course continue to satisfy the other requirements, like the Lorentz invariance 
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/ cluster separability requirement. 



Among these requirements, the absence of continuum dissolution is a technical necessity, if we want 
to get normalizable bound states. The correct nonrelativistic limit is of course a must. The possibility 
of improving the chosen approximation could enable us to estimate the precision of this approximation 
and possibly to satisfy the other requirements up to the considered order at least. 

Section 2 is devoted to the two-fermion problem, revisited in order to define the notations and present 
the building blocks and cluster-separated limits of our future three-fermion equations. In section 3 
we write directly a cluster separable 3D equation, by combining the relativistic free hamiltonians and 
the 3D potentials obtained in the reduction of the two-fermion Bethe-Salpeter equations, putting in 
each 3D potential the third fermion on the mass shell. Switching off two of the three mutual interaction 
potentials gives then a free Dirac equation for a spectator fermion plus the two-body equation which would 
be obtained directly by 3D-reducing the two-fermion Bethe-Salpeter equation, for the two interacting 
fermions. This insures the Lorentz invariance / cluster separability property. Positive-energy projectors 
solve the continuum dissolution problem. Our 3D equation can also be obtained from the three-fermion 
Bethe-Salpeter equation, by neglecting the three-body kernel and replacing the two-body kernels by 
positive-energy instantaneous (i.e. independent of the relative energies) approximations, equivalent at the 
cluster-separated limits. The remaining dependence (in the free propagators) of the approximated Bethe- 
Salpeter equation in the relative energies can then be eliminated by performing closed path integrals. This 
derivation of our 3D equation from an approximated Bethe-Salpeter equation enabled us to transform the 
neglected terms, known at the Bethe-Salpeter equation level, into correction terms to the 3D equation. 
We achieved this by performing three consecutive expansions of the correction terms: the three-fermion 
Bethe-Salpeter equation is first transformed into a set of three coupled equations for three wave functions 
depending each on one two-fermion total energy, then into a set of three 3D equations and finally into a 
single 3D equation. The resulting expression for the 3D equation is rather complicated, as it combines 
three series expansions, but we use it to write a manageable expression of the first-order corrections to 
the energy spectrum. The contributions of the residues of the poles of the propagators to this energy 
shift are easy to compute in terms of the Bethe-Salpeter kernels with two of the three fermions on their 
mass shell, in the spirit of Gross' spectator model Jl?], [l8|. Besides, we have also the contributions of 
the singularities of the kernels themselves to the closed path integrals. These contributions can not be 
a priori neglected. In section 4, we examine some other possible choices of the "approximation" of the 
free propagator, including Sazdjian's covariant propagator [^[ [l^, [l9|. In combination with a covariant 
substitute of the projector on the positive free-energy states, this propagator leads to a 3D potential given 
by a series in which each term is separately covariant, so that we can truncate the series without breaking 
the covariance. In section 5, we examine an alternative 3D reduction method built around instantaneous 
approximations of the Bethe-Salpeter kernels. The non-hermitian resulting 3D potential can be made 
hermitian either by performing a second series expansion, easily combinable with the series expansion 
generated by the 3D reduction, or by adaptating the starting approximation of the Bethe-Salpeter kernel 
(Phillips and Wallace's method [2Q] ). It is interesting to note that Phillips and Wallace's method and the 
propagator-based 3D reduction method give both the correct scattering operator (a bonus not implied by 
the starting homogeneous Bethe-Salpeter equation, designed for bound states), while the ordinary kernel- 
based 3D reduction does not. Section 6 is devoted to a comparison of our "hand-made" 3D equation 
with Gross' spectator model's equations |l7], [T^] in Faddeev's formalism. In this model, the dependence 
of the two-fermion transition operators in the relative energies is taken into account by putting different 
fermions on their positive-energy mass shells, according to the neighbours of each two-fermion transition 
operator in the expansion giving the three-fermion transition operator. Each of our two-fermion transition 
operators is everywhere approximated in the same way, so that our Faddeev equations are much simpler 
and can be transformated into our wave equation. The difference between both models corresponds to 
higher-order contributions, which we recover in our calculation of the higher-order terms, together with 
the supplementary contributions, not to be a priori neglected, of the singularities of the kernels to the 
closed path integrals in the relative energies complex planes. In section 7, we compare the one-high mass 
limits ("two-body limits") of our three-fermion equation with our previous two-fermion in an external 
potential equation [^6| . Section 8 is devoted to conclusions. 
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2 The two-fermion problem. 



2.1 Notations. 

We shall write the Bethe-Salpeter equation for the bound states of two fermions Q as 

$ = G°A"$, (1) 

where $ is the Bethe-Salpeter amplitude, function of the positions xi, x 2 or of the momenta pi,p 2 of the 
fermions, according to the representation chosen. The operator K is the Bethe-Salpeter kernel, given 
as a factor of the kernel of an integral equation in momentum space by the sum of the irreducible two- 
fermion Feynman graphs. The operator G° is the free propagator, given by the product G\G^ of the two 
individual fermion propagators 



G° = G?G°, G° = i fa = 2 P<0 \ hi ft (2) 



where the hi are the Dirac free hamiltonians 

hi = di .pi + ft m l (i = 1, 2). (3) 

The Bethe-Salpeter kernel K should contain charge renormalization and vacuum polarization graphs, 
while the propagators G° should contain self-energy terms (which can be transferred to K ^lj ]). 
In this work, we consider only the free fermion propagators in Gf and the "skeleton" graphs in K, and 
we hope that the inclusion of the various corrections would not change our conclusions. 
We shall define the total (or external, CM, global) and relative (or internal) variables: 

X = ~0i + x 2 ) , P = pi+p 2 , (4) 

x = xi-x 2 , P=-^{Pi-P2)- (5) 
and give a name to the corresponding combinations of the free hamiltonians: 

S = h 1 + h 2 , s = -{hi - h 2 ). (6) 
We know that, at the no-interaction limit, the system is described by a pair of free Dirac equations: 

(pio - fci)* = 0, ( P 2o - fta)* = 0, (7) 
where VP depends on x\,x 2 . Let us also write their iterated version 

(p? - E 2 )V = 0, ( P 2 20 - E 2 2 )^ = (8) 

with 

E i = ^=(f i +m 2 i )i. (9) 

Interesting combinations can be obtained from the sum and differences of the equations (|t]) or of the 
iterated equations (||): 

(P o -S)V = 0, (po-a)* = 0, (10) 

H°V = 0, (p o -M)* = (11) 

with 

H° = 2[{p\ - ml) + (p 2 2 - ml)] PQ=ll = P 2 - 2(E 2 + E 2 ) + 4 M 2 , (12) 

„ = -L(E? - El) = -L(ft? - A£) = £ (13) 

The choice (|]j5]) of the variables X and p was partially arbitrary: other linear combinations are possible, 
provided the canonical commutation relations are satisfied. In the nonrelativistic framework, the choice of 
X is weighted by the masses, and p is chosen accordingly. A frequent choice in the rclativistic framework 
is a weighting by the free energies, which leads to a simpler complementary equation po^ = instead 
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of (|Tl]), but has the drawback that the coefficients of the weighting are not constant. We think that our 
choice (|4]J^) leads to simpler expressions, even if taking the nonrelativistic limit is less straightforward. 

Our free equations ([lO]) are eigenvalue equations, the energies P and p being the eigenvalues and 
S and s the corresponding operators. Equation (0) with H° is a fourth-order equation in Pq, with 



four solution for a given set of values of the spatial momenta. Equation (11) with po is the eigenvalue 
equation of an operator \x, depending itself on Pq. In the Bethe-Salpeter equation ([!]), Pq remains a 
constant, but po becomes an operator in position space, or a variable in momentum space. The 3D 
reduction will consist in eliminating this variable by performing an expansion around a fixed combination 
of the other ones ( s or /.i in the following) . 

2.2 3D reduction of the two-fermion Bethe-Salpeter equation. 

In the zero-order approximation, the free propagator G° will be replaced by a carefully chosen expression 
G s , combining a constraint like 8{pq — ^) fixing the relative energy, and a global 3D propagator like 
— 2iir(Po — 5) _1 /?i/32- The argument of the S and the inverse of the propagator should be combinations 
of the operators used in the free equations (at last approximately and for the positive-energy solutions). 
There exists an infinity of possible combinations [8, 12, 20, 22-34]. The best choice depends on the 
quantities one wants to compute (energy of the lowest state, hyperfme splitting, recoil of a nucleus, etc..) 
and on the properties one wants to preserve exactly in the first approximation (cluster separability, 
Lorentz invariance, heavy mass limits, charge conjugation symmetry...). All choices would be equivalent 
if all correction terms could be computed, but this is of course impossible. 

We shall see below that the reduced wave function (from which the relative time-energy degree of 
freedom can be trivially eliminated) is given by 

= G" 5 (G°) _1 $. (14) 

The choice of the constraint fixes the 3D hypersurface (po = /i for example) on which we want to work. 
The remaining of G 5 is a purely 3D operator, the different choices of which result in different 3D operators 
applied on a common basic 3D wave function and in different rearrangements of a common reduction 
series giving the 3D potential. The various propagator-based 3D reductions of the literature can thus be 
classified according to the constraint they use. Once this constraint chosen, we can only write different 
equivalent forms (or sometimes projections) of the same 3D equation. It is the unavoidable truncation 
of the reduction series which makes the difference (numeric, if we simply want to compute predictions 
of field theory; more fundamental, if we want to write constraint theory equations with QCD inspired 
potentials). 

Two natural choices for the constraint are S(po— s), based on the first-order equations (flO|) and S(po—fi), 



based on the second-order equations (11). In S(j>q — s), s is an operator which can be diagonalized in 
momentum space using the four projectors on the subspaces corresponding to the different signs of hi, hi 
(see below). By contrast, fx is a momentum-depending number, which reduces to (mj- m^/ZPo in the 
two-fermion center of mass frame P = which can be defined in the pure two-fermion problem (but not 
in the two-fermion plus potential problem or in the three- fermion problem) . Other constraints could also 
be chosen, such as that of Gross ||, which puts one particle (normally the heaviest) on its mass shell. 



We made a nonexhaustive review in ref. 12 



If we consider the contributions of the poles of 

G° = -j -j Pifo (15) 

2-P0 + Po - m + iehi ^P - p a - h 2 + ieh 2 

in an expression like KG°K, we must perform an integration with respect to pq. If K is instantaneous, 
we get 

d P0 G°(p ) = -2inTg° Ptfo, g° = 1 . - (16) 

P - b + teP 



where 



can also be written 



T = \{ti + r 2 ), n = -A^ = Y- = ^s 11 ^)- 

A«=AjAj, Af = ^i. (18) 



G 



When K is not instantaneous, we must add the contributions of its singularities. Furthermore, in the 
residues of the poles of G° we must take K at pw = hi or at p 2 o = h 2 , according to the chosen integration 
path and to the sign of r. In order to preserve an explicit symmetry between the two fermions, we shall 
use the constraint S(pq — s), to be combined with g° : 

G ss (po) = -2nrTS( Po -s)g p 1 p 2 . (19) 

The operator r has a clear meaning in the basis built with the free solutions: it is +1 for hi, h% > 0, -1 
for hi, ho. < and zero when they have opposite signs. It comes from the dependence of the po integral 
on the signs of the iehi. 

The choice G s — G ss has two merits: a) It leads directly to Salpeter's equation without higher- 
order correction terms to the potential when the Bethe-Salpeter kernel is instantaneous, b) In the two- 
fermion plus potential problem and in the three-fermion problem, the operator r prevents the mixing of 
asymptotically free fermions with opposite energy signs, which is the origin of the continuum dissolution 
problem (see section 4 below). In order to reduce the number of terms at lowest order, we shall however 
keep only the A ++ part of r and choose 

G 5 (p ) = G sp ( Po ) = -2 i 7TA ++ S(po~ S )g°(3 1 p 2 . (20) 

We shall examine some other possible choices in section 4. 

Let us now write the free propagator as the sum of the zero-order propagator, plus a remainder: 

G° = G 5 + G R . (21) 

The Bethe-Salpeter equation becomes then the inhomogeneous equation 

$ = G°K^ = (G s + G R )K$> = * + G R K<$>, (22) 

with 

V = G 5 K<f> (= G S (G°)~ 1 $). (23) 

Solving (formally) the inhomogeneous equation (|2^) with respect to $ and putting the result into (p3|), 
we get 

^ = G s K(l-G R K)- 1 ^ = G S K T ^ (24) 

where 

K T = K(l - G R K)- 1 = K + KG R K + ... = (1 - KG R y x K (25) 

obeys 

K T = K + KG R K T = K + K T G R K. (26) 

The reduction series (Eq) re-introduces in fact the reducible graphs into the Bethe-Salpeter kernel, but 
with G° replaced by Cx. Equation ( p4| ) is a 3D equivalent of the Bethe-Salpeter equation. 
The relative energy dependence of eq. (|24|) can be easily eliminated: 

V = 5(p -s)^ (27) 

and ip obeys: 

$ = -2i7rg°A ++ J dp' dp S(p' Q ~ S )l3iP 2 K T (p' ,po)S(p -s)^. (28) 
Using the identity tp~A ++, ip, we can write 

i^=g°V^, V = ^Tr/^i^+^s), (29) 

(3il3 2 K T++ ( s ,s) = A++ / dp' dp 5(p' Q - s )pi/3 2 K T (p' ,p )6(p ~s)A ++ . (30) 



Note that we write (p' ,po) but (s,s), as we keep s in operator form. This operator can be diagonalized 
in the spatial momentum space by using the A lJ projectors. The eigenvalue will depend on the position 
of s in the formula: the eigenvalue of the first s in (|3^) will be built with the final momenta and that 
of the last s will be built with the initial momenta. 
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The inversion of the reduction is given by 



$ = {1-G R K)- 1 ^ = ( 1 + G R K T ) * = (l + G°K T -G S K T )V = G°K T ^ (31) 
or, explicitating the relative energy: 

*0/ ) = G a (p' )K T (p' a , s )^. (32) 

The splitting of G° into two terms containing a S is the origin of unphysical singularities in the 
terms of K T when the argument of the delta vanishes on the singularities of K. When the full K T is 
computed, the singularities of the different terms cancel mutually. When K T is truncated, some of the 
unphysical singularities have to be removed by hand |l7| . 

This 3D reduction can also be described in terms of transition operators. The 4D transition operator 

is 

T = K + KG K + • • • (33) 
and K T can be obtained by keeping only the G R part of G° in it. We have also 

T = Kil-^K)- 1 = K{l^G R K-G s K)- 1 = K{l-G R K)- 1 {l-G 5 K{l-G R K)- 1 )- 1 

= K T (l-G s K T )- 1 =K T + K T G 5 K T + ■■■ (34) 
so that the 3D transition operator 

T 3D = V + Vg°V + ■■■ (35) 

is also given by 

T 3D = -2iTT(3 1 f3 2 T ++ (s,s). (36) 

We see that the 3D transition operator is a constrained form of that of field theory. Both operators 
become equal to the physical scattering amplitude when both fermions are on their positive-energy mass 
shells. This was not guaranteed a priori, as our original two-fermion Bethe-Salpeter equation ([!]) was 
valid only for bound states. Our 3D equation (29) is a bound state equation too. To include the scattering 
states we should add an inhomogeneous term, or write the equation in the form 

(P -E 1 -E 2 )i> = Vil>. (37) 

Conversely: 

T = K T + K T G S K T + K T G S (K T + K T G 5 K T + ■ ■ -)G S K T 

= K T + K T G 5 K T + K T G S /3i/3 2 T 3D G S K T . (38) 

A Bethe-Salpeter equation leading directly to the same 3D reduction can be obtained by replacing the 
kernel K by the instantaneous positive-energy kernel K T++ (s, s). In this case the 4D transition operator 
T becomes equal to T++(s,s) = -f3 1 f3 2 T 3D /2iir. 

In contrast with the textbook cases, the potentials deduced from field theory arc in general energy- 
dependent (although one starts often with an energy-independent approximation). We shall thus get a 
set of eigenvalues W x (Pq) corresponding to mutually orthogonal eigenstates tp (Pq) of the hamiltonian. 
The physical energy spectrum will be given by the solutions of the algebric equations 

P = W x (P ) (39) 

and the corresponding wave functions are no more orthogonal with the the usual scalar product, as they 
are no more eigenstates of the same hamiltonian. The usual methods of perturbation calculation are to 
be revisited [ ^5| . 

The first-order correction term to the energy is 

P Q - P (0) = -2iir < A++ fcfa (KG R K)(s, s) A++ > . (40) 

bmce P (0) , (KG R K)(s,s) and the zero-order wave functions are in general total energy dependent, 
( f4(i| ) is in fact a numerical equation in P (or a matricial equation in case of degeneracy). If the Pq 
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dependence of the zero-order potential is of higher-order, one has only to take the Pq dependence of 
Pq into account. 

What is lost in our 3D equation compared with the starting Bcthe-Salpeter equation? The missing 
information is given by equation ([32]), which gives the bound state solutions of the Bethe-Salpeter equation 
in terms of the corresponding solutions of the 3D equation by adding the dependence in the relative 
energy. Furthermore, the transition operators are the same on the positive-energy mass shells, although 
the homogeneous Bcthe-Salpeter equation was originally deduced from the inhomogeneous Bethe-Salpeter 
equation (for the propagators) by postulating a bound state [[!]. Were some bound states lost in the 3D 
reduction, like the "abnormal solutions" related to excitations in the relative time degree of freedom? 
We have studied this still open problem elsewhere. Our conclusions were: 1). The abnormal solutions 
of the Bethe-Salpeter equation, when they exist, can also be obtained in the 3D equation via the total 
energy dependence of the potential, leading to multiple solutions of ( |39| ) for each value of A |3(| |37|. 2). 
The abnormal solutions found until now are artefacts due to the ladder approximation, and should not 
appear when all crossed graphs are included |37|, Q . 



3 The three-fermion problem. 

3.1 A cluster-separable three-fermion 3D equation. 

In the three-fermion framework, we shall append indices (ij)= (12) , (23) , (31) to the two-fermion symbols 
and write Py, ify, etc.... The index k=l,2,3 will denote the third fcrmion. Some of these symbols will 
also get a new three-fermion meaning when left unindexed, such as 

P = Pi+P2+P3, S = h 1 + h 2 + h 3 . (41) 
A cluster-separable three-fermion 3D equation can be obtained by combining three two-fermion potentials 

V> = 77^ ( V i2 + V 23 + V 31 )tJ> (42) 

given by our 3D reduction in the two-fermion framework. The potential V given by ( |29| ) depended on the 
total two-fermion energy Pq which we did not wrote as an argument, because it was conserved. When 
using these two-fermion potentials in the three-fermion framework, we shall replace the now nonconserved 
two-fermion energy argument Pjjo by the operator Pijo + (pko — hk) = Po — hk, which is the difference 
between the conserved three-body energy and the free energy of the third fermion. The potential Vij 
is thus no more strictly a two-fermion term, and the third fermion is no more strictly a spectator when 
there are interactions between the cluster (ij) and the fermion k. At the (ij)+k cluster-separated limit, 
however, the fermion k will be free and Pq — hk will be replaced in Vij by the now conserved P^o • Let 
us indeed "switch off' the potentials V 2 ° 3 and in equation (|42|). This equation can then be splitted 
into 

ip = ipi2ip3, Po = P120 + P30, (43) 
ijJi2 = -5 %, ^rVi2{Piw) ipu, (P30 - h 3 ) ^ 3 = 0, (44) 

-T120 — Hj\ — HjI 

i.e. equation ( p9| ) for the (12) subsystem and a free Dirac equation for fermion 3. Our 3D equation flUf ) 
satisfies thus clearly the cluster separability requirement. Furthermore, the three cluster-separated limits 
are exact equivalents of the corresponding two-fermion Bethe-Salpeter equations. This cluster separability 
is a property of the equation, or, more exactly, of the full Green function. For a given scattering solution 
it is also possible to take the cluster-separated limit at fixed P. This is not possible for the bound state 
solutions. 

In replacing Pyo by Po — h% we exploit the fact that the opposite cluster is at the cluster-separated 
limit a free fermion. With four fermions, we should introduce, for exemple, a potential V34 inside a 
potential V12 and vice-versa, getting an unmanageable infinite nesting of potentials. This limits our 
recipe to the case of three fermions. 

Our two and three-fermion equations were written until now in an unspecified reference frame. We 
shall consider three possible answers to the covariance challenge: 

- Write explicitly covariant equations. This is not the case ou our two and three-fermion 3D equations, 
built around non-covariant "approached" free propagators. We shall briefly present in section 4 a 3D 
reduction method based on Sazdjian's covariant "approached" free propagator || [l2| [l9|| . 
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- Choose the total rest frame P = as the reference frame. All quantities can then be made invariant 
by building invariants with the 4-vector P. 

- Search for an implicit covariance by equivalence with an explicitly covariant equation. A 3D equation 
obtained by reduction of an explicitly covariant Bethe-Salpeter equation will for example be implicitly 
covariant if all terms of the series generated by the reduction are kept, even if the terms of the series are 
not separately covariant. If the series is truncated, the 3D equation will be (implicitly) covariant up to 
a given order only. 

What about our equation (fl2|)? We can write it in the three-fermion reference frame. If we "switch 
off' the potentials V23 and V31, the equation for the subsystem (12) will still be written in the global 
three-fermion rest frame, but this link with fermion 3 will disappear if we go back to the original two- 
fermion Bethe-Salpeter equation, which is explicitly covariant. If the series generated by the reduction 
of this equation is truncated, we loose the exact covariance of the two-fermion equation and the exact 
cluster separability of the three-fermion equation. 

A well known difficulty in the writing of equations for two mutually interacting relativistic particles 
in an external potential and for three or more mutually interacting relativistic particles is the continuum 
dissolution problem: if we combine asymptotically free particles with opposite energy signs, it becomes 
possible to build a continuum of solutions at any given total energy. This renders the building of nor- 
malisable bound state solutions impossible. The introduction of our positive free-energy projectors is 
the usual way of solving this problem. Our 3D equation has solutions built entierely in the A +++ = +1 
subspace and satisfying 

^ = H IF— 7? F- ( Fl2 + v ™ + ( 45 ) 

JTq — £j\ — -E/2 — ft 3 

We shall come back to this point in section 4. 



3.2 Positive-energies instantaneous appproximation of the three-fermion Bethe- 
Salpeter equation. 

The Bethe-Salpeter equation for three fermions is 

$ = [G°G°A 12 + G° 2 G Q 3 K 23 + G%G%K 31 + G?G§G°A 123 ] $ (46) 

where A123 is the sum of the purely three-body irreducible contributions. Let us first consider the case 
of three two-body positive-energy, instantaneous and cluster-energy independent interactions Kf, (i.e. 
containing A^ + projectors and independent of p'ijo,Pijo, Pijo )• The Bethe-Salpeter equation becomes 
in this case 

$ = G?G°G° ft&ft [ V° V12 + V 2 % ^ 23 + V£ V 31 ] (47) 

Vg = -2m ^ijipm) = Pk Go, 1 J dp ij0 (48) 
Replacing <I> in ( f4S| ) by its expression (^) leads to a set of three coupled integral equations in the ipy : 

^12(^30) = ^ / d P™ G?G§ A/32 [ V& fefeso) + V° 3 ^(pio) + V 3 \ ^ifeo) ] (49) 
where pw,p 2 o must be written in terms of Po,P30,Pi20 : 

P10 = 2 h Pl20 ' P20 = 2 Pl20 ' ^ ' 

and similarly for ip 23 and "031. We shall now search for solutions ipijiPko) analytical in the lower 
Im(pfco) <0 half plane and perform the integration ( [l9"| ) by closing the integration paths around these 
regions. The only singularities will then be the poles of the free propagators. The result is 

V>i 2 (P3o) = 775 ~ [V? 2 Mpw) + V&M^i) + V 3 °^ 31 (E 2 ) ] (51) 

(r - £>) - [p 30 - n 3 ) + le 

and similarly for ip 23 and ip 3 i (from now on, we shall often write the equations for the (12) pair only, 
understanding the two other equations obtained by circular permutations on the ijk indexes). Solving 
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( |5l"| ) with respect to "012 (P30) gives a wave function which is analytical in the Im(p 30 ) < half plane, and 
confirms thus the existence of such solutions. Furthermore, equation (^TJ) shows that the three projections 
A^ipij(Ek) are equal (let us call them ip) and satisfy the 3D equation 

^ = p JF [ ^1° + V & + < ] 1>- (52) 

Fq — Oil — Ui-2 — HiZ 

Moreover, t/j is the integral of the Bethe-Salpeter amplitude with respect to the relative times: 

?A = 7r-~ (dpo® = — [ dp 10 dp2odp 3Q 5(p w + p 2Q + p 3 o - Po) (53) 

Let us now choose 

K% = K^++( Sij , Sij ,P -h k ) (54) 



The 3D equation (|5^ ) is then built with the two-body potentials V,^ = Va (Pp—hk) , with the Vij defined 
by ( ^9|) in section 2 for each fermion pair, and becomes equation (|45"1). 

Let us now consider a non-approximated three-fermion Bethe-Salpeter equation (|46| ) and try to trans- 
form it into our 3D equation ( |52|) above, completed by higher-order correction terms. It first seemed us 
a good idea to start with approximations giving the 3D equation (js^) in a more direct way than above, 
like 

G?G° « G?f (55) 

G?G° <5(p 3 o - P30 ) « Gf^ «5(p 30 - /i 3 ) (56) 

Ku S( P ' 30 - p 3 o) w A 3 + ° 2 <5(^ - S 3 ). (57) 

Each of these approximations gives in each case the 3D equation (f52j), but the correction terms spoil 
the result: with ([55|) or (|5^) they contain supplementary zero-order contributions. With ([57]) they 
contain no supplementary zero-order contributions, but the correction terms can not easily be sorted by 
increasing order. We must thus proceed more carefully. We finally decided to start with positive-energy, 
instantaneous and cluster-energy independent approximations of the two-body kernels, like ( |54| ) or its 
Born approximation, adding the third A^~ projector in order to avoid unnecessary complications, and to 
proceed in three steps, inspired by the demonstration (|4^)-([52|): 

i) Transform the Bethe-Salpeter equation (^6|) into a set of three coupled equations for the ipij(pka)- 

ii) Transform this set into a set of three coupled 3D equations for the ipij(Ek) (these are not the 
Faddeev equations). 

iii) Transform this last set into a single 3D equation for ip, defined as the mean of the ipij(Ek)- 

3.3 Three coupled equations for three cluster-energy depending amplitudes. 

Let us now consider the unapproximated Bethe-Salpeter equation ( |46| ) and define 

G° = G? G° 2 G° 3 , (58) 

K = (G*)- 1 ^ + (G°)~ 1 K 2 3 + {G 2 ')~ 1 K 31 + K 123 . (59) 
The Bethe-Salpeter equation ([l6]) takes then the compact form: 

$ = G°K$. (60) 

Let us also define 

K° = (Gg)- 1 ^ + (G?)- 1 ^ + {G 2 )~ 1 K° 1 , (61) 

K R = K - K° , (62) 

each Ky being independent of all relative energies and containing A +++ projectors (in order to render 
the following calculations simpler, we now add the A^ projectors to the A+ + projectors of the previous 
section). Let us now perform the transformations 

$ = G K $ + G°K R <£ (63) 

*= {1-G K R )- 1 G°K°<P = G K K°$ (64) 
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where 

G K = G° + G°K R G° + G°K R G°K R G° + ■ ■ ■ = G° + G KR . (65) 

Since eq. ( |64| ) ends with K°&, we can transform it into a set of three coupled equations for the ipij(pko), 
now defined with the A +++ projector: 



We get 



(Pko) = A+++ fa Gg fc J d Pij0 $. (66) 

^12(^30) = X! / d P*oGl2,»j(P30>PM)-K$ ipijiPko) (67) 
Gf 2>i >30,Pfco) = A+++ ft [G° 3 {p' 30 )] - 1 [ dp' 120 dp ij0 G K (p' , Po ) (3 k (68) 



with 



3.4 Three coupled 3D equations. 

Let us now compute the contribution of G° to (|68|). We get (omitting the A +++ already contained in 

K°y. 

G? 2 ,i 2 (P3o>P3o) = 775 e\ — T~t \ i ■ fafaHPso -p 30 ) (69) 

Gi2,a3(Pao,Pio) = ft G?(Pio) ^(Pq -Pio-p 3 o)ft (70) 

G? 2 , 3 l(P30 5 f20) = A G°(p 20 ) G°(P - P20 - P 30 ) A (71) 

The product of propagators ( |70| ) will be included in a clockwise path integral in p±o. With positive-energy 
instantaneous kernels, this integral is given by the residue of the pole of G5(pio) at piq — E\. This is no 
more true in the general case, but we can still isolate a part of G?(pro) which leads to the same result: 



with 



so that 



G°(p 10 ) = Gf(p 10 ) + Gf (p 10 ) (72) 

Gf(p 10 ) = -2<7r<S(p 10 -f7i)/3i (73) 
Gf (pio) = ^ ft (74) 



(P -5) - (p 30 - P 3 ) +«e 



and similably 
Putting all together, we have finally 



— 2iir 

<?i2,3l(P30>P20) = 773 m T7 F N T ■ ftft £(P20 - # 2 )- (76) 



— 2l7T 

= (P -S)-(pi 30 -E 3 )+*e ^ ^^(P'so) 

— 2l7T 

[ft/% # 2 3 fe(Pl) + ftft if 3 °i ^3l(P 2 ) . 



(P -5)-(p 30 -£ 3 )+ie 

+ J d Pl0 fa G R ( Pl0 ) G°(P - p 10 - p 30 ) A K 2 ° 3 febio) 

+ / dpm fa G R ( P20 ) G?(P - p 20 - p 30 ) ft ^ 3 °i ^si(p 2 o) 

+ £ / **o Gf 2 ^ (p 30! Pfe0 ) K° ^ {Pko) (77) 

If we neglect the last term of ( |77j ) the two preceding terms vanish also by integration, and we recover the 
result of the positive-energy instantaneous approximation. We shall thus move the first term of (|77|) to 



12 



the left-hand side, consider the two next terms (second line) as the principal contributions, and the three 
last terms as perturbations. The difference between this approach and that based on ( |57| ) lies in the fact 
that we do not replace immediately p' 30 by E3 in the first term. We get 

Mp' 30 ) = 9i2(p' 30 ) [t^lM^i) + V&iPsiiEz)] 

+ 9M0) l9°Mo) r 1 E / G i2 T Mo,Pko) K% ^-(pfco) (78) 

k 

where Ofay is defined from G^lj by including the 3th and the 4th lines of (|77]) in the last one, and 



3i2 (P30) = 75 o — t7o \ 1 1 1 em 1 • ( 79 ) 

[P Q -S- V&) - (p 30 - E 3 ) + it 

(P - S) - (p 30 - E a ) + ie 



The iterations of ( |78| ) lead to 

^12(P') = 512 (P30) [^23 (Pi) + ^31 (P^ 



2 (p 3 o ) [ 9i2 (Pao ) ] 1 E / dpk0 G ^ T zj (P30 > PkO ) Kfj 



912 

k 

9iMo) [ V? k ^k(-B<) + V$ Vfci(^) ] • (81) 
defining G TT as the sum of the iterations of G KT : 

GfZMoiPko) = Gi^04).P*o) + E / dpUGEh'^oMK^ 

k> 

ft'i'(P*'o) [^(ffc'o)]" 1 ^^^'^,^) + • • • (82) 

We can now take (^lj) at p'=p3, multiply at left by [ <?i2(p3)] _1 , and move the term in ^2^12(^3) 
back to the right-hand side: 

^12(E a ) = 1 —— [V^12(E 3 ) + K>° 3 V23(Pl) + V^ 31 {E 2 )} 

+ E M ^ [ V Jk 1>ih ( E i) + V ki V>« (Ej ) ] (83) 

with 

M 12 ,ij = J dp k0 G? 2 T tJ {E 3:Pk0 )K° 3 -gijipko). (84) 
3.5 One 3D equation. 



Let us write (|33|) and the two other equations obtained by circular permutations in matricial form. We 
shall define 

* fc = M E k), n« = i n = l-n. (85) 

The matrix II is the projector (II 2 =11) on the u =1 vector (of norm 3). We have 

= tpu, if; = - (* x + * 2 + * 3 ). (86) 
3 



We can write eq. (B3J) in the form 



* = Mtf, M = M°+M , M oa = Vo ,, etc... (87) 

Po - S + ie 
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The action of the projector II on M° is 



II M° = M°, M° n = — V°, V = V.% + VS, + K , . (88) 

P - S + 16 12 2i JL \ I 



We shall now transform the matricial equation (|87j) into a scalar equation for ip. The expansion of ([87]) 
with respect to IW gives 

n * = um(i + um + Tim Tim + ■ ■ •) n * (89) 

in which ILM =JIM R \ Taking the scalar product with u leads finally to 



1> 



' V° + \u T ( M R + M [ UM R + I1M R I1M R + ■■■}) u 



P -S + ie 3 



tp. (90) 



In this equation, V° contains the projector A +++ at left and at right, while M R contains it at right. 
We can thus trivially transform (^0|) into an equation for A +++ ip. 

3.6 First-order energy shift. 

The 3D reduction above is rather complicated, since it implies three levels of series expansions. The 
computation of the first-order energy shift remains however tractable. The energy shift is, keeping only 
the terms which contribute at first-order: 

P - P (0) = < ( p o - S) \u T ( 1 + M°TT ) M R u > . (91) 

We can make the replacements 

1 + M°Tl = 1 + M - 11 V° -> M° (92) 
as the contributions of the first and third terms cancel mutually. We remain with 



Po - P (0) = < (P - S) - u T M° M R u > 



= J2<V$ jf M i , j ,, ij (V -V? j )> (93) 

k',k 

The contribution of the terms with i'f = 12 is then 

A12 = < V° 2 I d PM (E 3 , Pk0 ) K% m ( Pk0 ) (V°-V%)> (94) 
k J 

At first order, we can replace Gj^ij by the non-iterated Gy^ij an( l further by Gf 2 %, as the 3th and 
4th lines of ( \n\) will lead to null integrals. Using the definitions (^3j|), ([58]) we get 

A 12 = < V? 2 ]T I %30 - E 3 )dp' d Po fa [GUp'^V 1 

k J 

G KR (p' ,Po) Pk K% 9ij ( Pk0 ) [V-V$)> (95) 

Replacing V by P -S and G KR by G°K R G° leads to 

Po - P (0) = ^ < fdp' d Po V F {p' ) G (p' )K R (p , ,p )G°(p Q ) fcfcfc VM > (96) 

with 

Vi(po) = V °o 9ij (Pko) g-/(E k ) (97) 



14 



V F (p' ) = -2^ V h> <*(pU ~ K>) fa [Gl> (P'k'o)} l - (98) 
fc' 

In writing (|98|), we mean that [G°, (p^q)] 1 cancels the G k , {p' k , ) contained in G°(p' ) before the 
replacement of p' k , Q by E' k ,. The energy shift ( p6[ ) contains 36 terms: there are 4 terms in K R (if we 
do not forget the irreducible three-body term) times 3 terms in Vj and 3 terms in Vp . Our reason 
for introducing these Vj and Vp is the fact that they both could be replaced by V° if K R were 
instantaneous. 

The singularities to consider when performing the integrations are the poles of the propagators, the 
singularities of K R and the poles of Vj. It is always possible to choose the integration paths, in each term, 
in order to avoid contributions from the poles of Vj. The contributions of the poles of the propagators 
to the energy shift @) are then easy to compute, as gij(pko) in ( p7| ) is to be taken at pko = hk- We get 

R12 =<E v$i> v if a v "j > (") 
fc'fc r ° b n b 

plus the similar R23 and R31 and a contribution i?i23 of the irreducible three-body kernel. By V^2 ba 
we denote — lit: (5\fiiK R 2 with the initial fermion a, the final fermion b and the spectator fermion 3 on 
their positive-energy mass shells. We have to consider 9 terms in 4 different groups: 



fc' = fc = 3 

fc' = lor 2, fc = 3 

fc' = 3, fc = lor 2 

fc' = 1 or 2 , k = 1 or 2 



6a = 11, 12, 21 or 22 
ba = fc'l or fc'2 
ba = 1 fc or 2 fc 
6a = fc'fc 



We still have some freedom in the choice of 6a which we could use to get a real energy shift, for example 
by choosing a = b in the three first groups. The contributions of an instantaneous part K R 2 of K R 2 
would sum into a < > but we could also include this K R 2 directly in K\ 2 . The corrections 

corresponding to the differences between the various V^ ba are a step towards Gross' spectator model 
(see below, section 6 and 0, p^|). 

The principal task remains however the calculation of the contributions of the singularities of the 
K R , for the integration paths defined by our choices of a6 above. As we shall see in section 5.3, these 
contributions could be as important as these coming directly from the relative time dependence of these 
kernels (they could however be partially cancelled by the contributions of two-body non-ladder terms 
and three-body terms ^ [l2[ [lTj ). A general expression of the contributions of these singularities can be 
obtained by splitting the propagators containing the poles considered in (B9J) , using (|72|), the all-G" 5 part 



giving (99) 



For more explicit results, we should specialize to specific three-fermion problems. It must be noted 
that the ratio gij(pko) g^i^k) of ( |97| ) is no more 1 in the terms containing G R (pko)- If the singularities of 
K R in the complex pko plane lie far from the positive-energy mass shell, this ratio can be approximated 
by —{pkQ—Ek)^ 1 g^{Ek) and takes a relatively low value. This suggests us to neglect the contributions 
of the terms containing G R (pko) to (p6|), as these contributions contain at least one large denominator 
more than the contributions of the terms containing G s k (pko) ■ This leads to the replacement 

Vi(po) - = -2tn J2 V?j Pk [GKpw)] ^S( P ko - E k ) + ]T V% (100) 

k^u k—u 

in (Q) . We isolated the contributions of the unconnected terms fc = u (like fc = 3 combined with fc' = 3 
and K R 2 )■ The resulting energy shift is now manifestly real. 



3.7 Cluster separability. 

Physically, if we "switch off" the (23), (31) and (123) interactions, we remain with a pair of interacting 
fermions (12) and a free fermion (3). We would like to read this on our final 3D equation, however: 

- The cluster-separability condition may be satisfied by an exact 3D reduction, but not by its truncated 
forms. 

- Even with an exact 3D reduction, the cluster separability property is a welcome bonus, not a 
requirement: the Bethe-Salpeter equation ([l6"|), from which we started, is indeed correct for three-fermion 
bound states only. 
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In our calculations above, the cluster-separability property is in general spoiled by our iteration of 
(|78|), which starts with a term which vanishes at the (12) cluster-separated limit. The (12)-limit equation 
is then built with the approached —21^ (5\j5 2 K\ 2 only. The cluster separability can however be restored 
by performing a supplementary iteration, like our replacement of V° by Pq — S in ( 93ft , leading to 
the cluster-separable correction (^9|) . Another possibility consists in our previous choice ( 5J) of the . 
With this choice, the cluster-separated limits are exact. Moreover, the corresponding two-body scattering 
amplitudes are also correct (of course, we do not mean that we should keep all the terms of the 
in an actual computation: we should only keep the ones which contribute up to the desired order). 
Simpler choices (as the Born approximations) would however not be uncorrect in the calculation of the 
three-fermion bound state spectrum. 



4 Other propagator-based 3D reductions. 
4.1 Two fermions. 



We have already seen two different choices of G s : Salpeter's propagator G ss (19) and the positive- 
energy propag ator G 5P d20l), which we used through this work until now. We could also use Breit's 
propagator fl39| : 

G 5B (p ) = -2m6(p -s)g°l3 1 /3 2 . (101) 

in the two-fermion problem only, as it would lead to continuum dissolution in the two-fermion in an 
external potential and in the three-fermion problems. The first-order constraint S(po~s) could also be 
replaced by the second-order constraint d(po—[i), which leads to a simpler one-body limit, especially with 
Breit's propagator (see section 7 below). In the spirit of Gross'spectator model, we could also choose the 
unsymmetrical or symmetrized propagators 

G 5G1 { Po ) = -2iirA+5{p2o-E 2 )g°p 1 p 2 . (102) 

G SGS ( Po ) = -2i7r| [AtS(p 20 -E 2 ) + A+Sip^-Et)] g° frfo. (103) 

in the two-fermion problem (for the three-fermion problem, see below section 6). 
With Sazdjian's covariant propagator H, |lj, |l^] 



G sz = -2in bi -7i +mi)( P2 . 7 , + m 2 ) 
Pi + V 2 - \ m i + m 2 ) + i e 



= -2nr^^£(po-M)<7°A/?2. (104) 

the 3D potential is given by a series of separately covariant terms, so that we can truncate it without 
breaking the covariance. In order to preserve this property in the two-fermion in an external potential 
problem and in the three-fermion problem, without meeting continuum dissolution, we must however 
combine this propagator with a covariant substitute of the noncovariant A ++ . A possible choice is the 
product 

e ++ = d(p 1Q )9(p 2 1 )d(p 2 o)d(p 2 2 ) = d(p 10 -\p 1 \)9(p 2 o-\p 2 \) (105) 

in which the constraint will replace pio by ^Po i /i. The sign of pio is indeed invariant when p 2 > 0. 
Another difficulty comes from the (Po + 5)/2|Po| operator which obliges us to perform supplementary 
manipulations after the 3D reduction in order to render the potential hermitian, so that the Born term 
itself is finally given by a series. We examined elsewhere [|40"[ how to use Sazdjian's propagator in the 
two-fermion, two fermions in an external potential and three-fermion problems, leaving the reader decide 
wether the covariance of the approximations balances the supplementary complications. 
With Salpeter's propagator, the 3D potential will be 

V s = -2iTTTp 1 p 2 K TS (s,s)T 2 (106) 

as T 2 ip S —ip S > with K TS built with G ss like K TP with G sp . This potential is hermitian in the r 2 = l 
subspace at fixed P with the scalar product 

(^ , ) = / d 3 p (p) r(p) (p) . (107) 
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4.2 Three fermions and continuum dissolution. 



The results of a 3D reduction of the two-fermion Bcthc-Salpeter equation can be used in the three-fermion 
problem at two levels: by choosing the corresponding approximations of the two-fermion kernels in the 
Bethe-Salpeter equation ( |46| ) , or by combining directly the corresponding potentials into a 3D equation 
(fl2|). Both approaches turn out to be equivalent in the case of our 3D reduction based on G sp . With the 
other choices, the 3D equation obtained by reducing the approximated Bethe-Salpeter equation would 
contain supplementary terms. In all cases, we get an exact three-fermion equation by "switching off' the 
interactions with the third fermion. However, the three-fermion 3D equation built with the potentials 
V B or V sz , based on approximated propagators without projectors, suffers of continuum dissolution. 
It is indeed possible to build a continuum of solutions with any a priori given total energy by combining 
asymptotically free fermions with opposite energy signs. Any physical bound state is thus mixed with such 
a continuum and the building of normalisable bound state wave functions becomes impossible [13-16]. In 
the pure two-body case the energy of a system (+, — ) in the total rest frame is 



E 1 -E 2 = JfP + m - JjP + ml = 7=rr=5 (108) 

and lies thus between mi — 7772 (whichever the sign) and zero. We have thus no problem if we make the 
assumption that the energies of the bound states lie between | m\ — m 2 | and mi + m 2 ■ Below mi — m 2 \ 
we meet the "strong field" problem. The strong field problem for the one-body plus potential and two- 
body systems, the continuum dissolution problem for the two-body plus potential and three-body systems 
are both consequences of the possibility of pair creation. 

In the three-body case, the mixing with asymptotically separated (12) (3) subsystems of opposite 
energy signs is excluded by the total momentum conservation, as in the two-fermion case. For the mixing 
with three-fermion asymptotically free states, let us consider for example 



E 3 - E x - E 2 = ^0* + m\ - ^{ +m{- ^Jf 2 + mj , p x + p 2 + p 3 = 0. (109) 

There is no lowest value (we can for example have p 3 = and pi = — p 2 arbitrarily large) . The highest 
values are obtained when p\ and p 2 have the same direction, and for pi/m\=p 2 /m 2 = — p 3 /(m\+m 2 ). 
In this case we have 



E 3 -E 1 -E 2 = ^f 3 + m 2 3 -\JpI + (mi + m 2 ) 2 (110) 

so that E 3 — Ei — E 2 lies finally between —00 and (7773 —mi —7772)0(7713 — mi —m 2 ). Symmetrically, 
Ei+E 2 — Es lies between (7771+7712 — 7713)0(7713— mi— m 2 ) and 00. If we assume that the energies of the 
bound states lie between (7771+7712-1-7773— 2Inf(mi)) (i.e. above the highest negative-energy threshold) 
and (mi+7772+7773), there is no degenerescence with the "one plus-two minus" states. On the contrary, 
no weak field assumption could prevent the degenerescence with the "two plus-one minus" states. 

The continuum dissolution problem is avoided with the potentials containing the positive-energy 
projectors A^ + , or their covariant substitutes . The 3D equation built with Salpeter's potentials 

, which include a Tjj at left and a rfj at right, can be splitted into two uncoupled equations in the 
TlT2T3 = ±1 subspaces and the solutions are eigenstates of tit 2 ts and free of continuum dissolution. The 
TlT2T3 = +1 subspace includes the three-plus and the two-minus one-plus components. The interaction 



term is hermitian for a scalar product built with T123 = sign (ti +t 2 +T3). Equation (42) becomes then a 
logical three-fermion extension of Salpeter's equation and exhibits a particle-antiparticle symmetry. We 
do not expect however to get a large contribution from the two-minus one-plus components. 



5 Kernel-based 3D reductions. 
5.1 Two fermions. 

Our 3D reduction of the two-fermion Bethe-Salpeter equation was built around an approximation of the 
free propagator. However, the first step of our 3D reduction of the three-fermion Bethe-Salpeter equation 
is built around instantaneous approximations of the kernels, while the second step is again built around 
approximations of the free propagators. It is therefore interesting to examine also how a 3D reduction 
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of the two-fermion Bethe-Salpeter equation can be built around an instantaneous approximation of the 
kernel. Let us thus consider an approximation K z of the Bethe-Salpeter kernel: 

K = K z + K R . (Ill) 

The Bethe-Salpeter equation becomes 

$ = G Q K Z <& + G°K R <P (112) 

$ = (1 - G°K R )~ 1 G°K Z 3> = G K K Z <$> (113) 

with 

G K = G° + G°K R G° + G° K R G° K R G° + ■■■= G° + G KR . (114) 

If we now specialize K z to an instantaneous positive-energy kernel (K z independent of po and equal 
to /3 1 /3 2 A ++ /3 1 /3 2 A" Z A + + ), eq. ( p!3|) leads to the 3D equation 

$ = {g Q + g KR )V z <j ) (115) 

with 

= A++ J dp $(po), V z = -2nr li x li 2 K z , (116) 

« 7 *« = -J-A++ fdp' dp G KR (p' ,p )l3 1 (3 2 A++. (117) 

2l7T J 

The first-order energy shift is, with the replacement g KR V Z (j> « g KR (Po — S)cf> : 

P - P (0) = ~ < (P°S) f d P ' dpoG (p' Q )K R (p' 07 p )G°(p )(3 1 (3 2 (P°~S) > . (118) 



Here again, all elements of ( |118j ) are tota l-energy dependent. Although (118) is symmetric, the pertur- 
bation potential (Pq~ S)g KR V z of (115) is not. This does not exclude the possibility of a real energy 
spectrum: if we admit an energy dependence of the potential (which we cannot avoid), we have also 
the possibility of performing an infinity of rearrangements of the equation. We can for example render 
the potential symmetric by treating g KR as a perturbation of g° (like G R as a perturbation of G s in 
section (2.2)). We can equivalently start with the corresponding 4D equation (113): 

$ = (G° +G KR )K Z $ = (f> + G KR K z $ (119) 

4> = G°K Z <£ = G°K z (l - G KR K Z Y X 4> = G°K K (j) (120) 

where 

k k =K z + K z G KB. R z + , . . (• 121 - ) 

is a symmetric positive-energy instantaneous potential, since it begins and ends with K z . The 3D 
equation is then 

V = 9°V K V, V K = -2zn 13x132 K K , V =[dp ^(p ) (122) 



and the first-order energy shift is still given by ( 118 ). As for the propagator-based reduction, we can 
represent the various contributions by Feynman graphs. We can indeed get K K by considering the 
transition operator 

T = {K z +K R ) + {K z + K R ) G° (K z + K R ) + ■■■ (123) 

and skipping some terms according to the two rules: 

- Each term must begin and end with K z 

- There must be at least one K R between two consecutive K z . 
The 3D transition operator is now —2iir(3if32T K with 

T K = K K + K K G°K K + ■■■ = K K (1 - G°K K y 1 

= K Z {1 - G KR K z )-\l - G°K Z (1 - G KR K Z )~ 1 )~ 1 
= K z (l - G KR K Z - G°K Z )^ 1 = K z (l - G K K Z )- 1 
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= K z (l-(l-G a K R r 1 G°K z )- 1 
= K z {\ - G°K R - G°K Z - G°K R ) = K z (l - G°K - G°K R ) 
K Z (1 + G°T)(1-G°K R ) 
Replacing then K R by K-K z and T(1-G°K) by K gives 

T K = K z + K Z G°K Z + K Z G°TG°K Z . 



(124) 



(125) 



In contrast with what happens in the propagator-based reductions, the on mass shell restriction of this 
operator is not the physical scattering amplitude of field theory. In fact, the Bcthe-Salpeter equation (|l|), 
without inhomogeneous term, is a bound state equation, and for the bound states we are only interested 
in the poles of (|l25| ). These poles are the poles of T, unless the integrations on p' ,po make some 
residues vanish. Getting also the correct physical scattering amplitudes was a welcome bonus of the 
propagator-based reductions. 

If we want T in terms of T , we can write (124) as 

T^l-GV)- 1 = K Z (1 + G°T) = (K -K R )(1+G°T) 

= T — K R (1 + G°T ) = (1 - K R G° )T - K R 
T = (1- K R G Q )~ l T K (l-G ^)- 1 + (l-A'^G )- 1 
= (1+T R G°)T K (1 + G°T R ) 

with 



K 1 



T 



T R = K R + K R G°K R + ■ 
We have a large freedom in the choice of K z . We could choose 



K 



/3 1 /3 2 A ++ /3 1 /3 2 X( S , S )A +H 



(126) 

(127) 
(128) 

(129) 



We could also keep only the ladder term of K, or a part of it in (129), like the Coulomb part of the 
one-photon exchange contribution. We could use K T (in this case we expect that the remainder of K K 
does not finally contribute to the total energy spectrum) or an approximation of it (one and two photon 
exchange terms, for example). Phillips and Wallace |2(j] suggested to choose K z in order to annihilate 
g KR exac tly or up to a given order. If we choose to annihilate it at first-order, the first-order energy 
shift ( 118 ) vanishes and we get 



K z = 



P 1 p 2 (g°)- 1 ^ ++ / dp' dpoG (p' )K(p' ,p )G (p )p 1 p 2 (g°)- 1 A ++ . 



(130) 



The skipping of terms in (123) is then improved, as we must now have at last two K R between two 

consecutive K z . If we choose to annihilate g KR exactly, we get, using (127), (117) and the relation 
qKR_ qo^Rq^ ■ 



-1 



(27T)* 



ftft ( fl F°)- 1 A++ / dp / dp G (p' )T{p' 0l p Q )G (p )p 1 p 2 ( ff )- 1 A + + = T 



K 



(131) 



which is the A ++ projection of the expression given in [Q (they use Salpeter's propagator and write 
thus A ++ — A instead of A ++ ). In this case, T K gives the correct physical scattering amplitudes 
when P° — >Ei+E2, as the operators (g )^ 1 kill the contributions of the singularities other than the 
poles of G°. 



5.2 Three fermions. 

We could also use the results of the kernel-based 3D reduction of the two-fermion Bethe-Salpeter equation 
as the starting point of a 3D reduction of the three-fermion Bethe-Salpeter equation, by choosing 

K% = K* (s,s,P -h k ) (132) 

as starting approximations in (|54|), in order to recover, at the 2+1 separated clusters limits, the kernel- 
based 3D equations. 
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5.3 One-photon exchange. 

Let us now compute ( |130[ ) with K given by the one-photon exchange graph, in Feynman gauge: 

_ 2ie 2 7i-7 2 
K ~ (2^)3 (p>- p )2+ie 



We have 



With 



rZ = ~ 7^4 A++ 7i ■ 72 A++ / (134) 



(2tt) 5 



J = J dp' w dp w 



P -E[- E' 2 

-E' 2 + ie' 2 ) 

En 

(135) 



(pi - £?i + zei) (Po - Ko - ^2 + ^a) 
1 P -E 1 -E 2 



(p'io - Pw) 2 - k 2 + it] (Pia - -Ei + *ei) (-Po - Pio - J5 2 + ie 2 ) 

and k=p'—p. Closing the integration paths clockwise, and taking the t\ — ■> limit before the ?7 — > 
limit in all terms, we get 

I = -4tt 2 (1 + Jl-Rk), (136) 

[E[ - Ei) 2 - k 2 2\k\ 

R = (E[ - E 1 - \k\)(P - E 1 - E 2 ) _ (P(-g 1 + |fc|)(P -^-£ 2 ) 

Pa - E[ - E 2 - \k\ P -E 1 -E' 2 -\k\ 

The first term comes from the residues of the poles of the propagators only, while the terms of Rk 
combine the residues of one pole of the propagators and one pole of the kernel. Modifying the integration 
paths or the order of the limits would give the same result (after rearrangements). At lowest order, we 
get 

R K « 2P a -E[-E! 2 -E 1 -E 2 . (138) 

The contribution of Rk to the bound state energy will finally be in a 3 , while the replacement of 
(E[ — Ei) 2 — k 2 by simply — k 2 (Coulomb potential) would lead to a correction in a 4 . This means 
that the relative energy dependence of the kernel, at the ladder approximation, contributes mostly by 
providing supplementary poles inside the integration paths. Beyond the ladder approximation however, 
the contribution of the first crossed graph will cancel the leading term in a 3 , leaving a contribution in 
a 4 ^2|. This is not a surprise, as it is well known that, in well chosen propagator-based reductions, all 
higher-order contributions cancel mutually at the one-body limit [pi n2l . 



By rearranging ( 136 , 137 ), we can also write 
. -4^ 2 / 1 



2\k\ \P -E' 1 -E 2 -\k\ Po-Et-E^-lk] 



(139) 



which is the expression obtained in time-ordered perturbation theory [^l| . Here, each term describes the 
emission of a photon by one fermion, followed by its absorption by the other one. 

5.4 Time-ordered perturbation theory 

In time-ordered perturbation theory (TOPT), one takes the Fourier transform of the full propagator (as 
given by the Feynman graphs method) with respect to the energy variables. Each graph becomes then 
splitted into several time-ordered graphs, and one integrates the result with respect to the time intervals. 
This leads to TOPT "Feynman rules". The TOPT graphs can be arranged in order to get an iterative 
equation for the propagator, and, by isolating the residues of the three-body bound state poles, to get 
an homogeneous 3D equation for the three-body bound state wave functions (this imitates the usual 
derivation of the inhomogeneous and homogeneous Bcthe-Salpeter equations). The cluster-separability 
property is translated by the factorization of the propagator at the different vanishing interactions limits. 
The two and three-body 3D potentials are directly given, at any order, by the TOPT "Feynman rules". 
Up to now, this approach has been applied to a very simple model: the exchange of a scalar meson 
between two or three positive-energy bosons or fermions at the ladder approximation |4l| . Our results in 
this work can be applied to more general interactions but should lead to more complicated calculations 
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than TOPT at higher-orders. A comparison of our first-order contributions computed with the simple 
Bcthe-Salpeter kernel used in j4l[] should be interesting. 

It must be noted that our motivation for introducing the positive-energy projectors is not that of 
[41|: we make the choice of approaching the exact equation for the positive-energy components, while 
[41 1 simply neglects these components as a temporary simplification. It will be interesting to see in which 
form the continuum dissolution problem appears in TOPT when the negative-energy components are 
taken into account. 



6 Relation with Faddeev formalism and Gross' spectator model. 

The nonrelativistic Faddeev equations can be obtained by transforming Schrodinger's equation (in a first 
step, without three-body terms) into a set of three coupled equations for three parts of the transition 
operator ( denotes the contribution of all graphs beginning by a (ij) interaction). The input is the set of 
the three two-body transition operators and the resulting series expansion contains only connected graphs 
(never twice the same two-body transition operator). This formalism is well adapted to the description 
of the various scattering processes, such as (12) + 3 — > (12) + 3 (elastic scattering), (12) + 3 — > (12)* + 3 
(excitation), (12) + 3^1 + (23) (rearrangement), (12) + 3 -> 1 + 2 + 3 (breakup). The (123) bound 
states correspond to the poles of this transition operator. 

The structure of these equations can be generalized to relativistic equations which are not necessarily 
(exactly) reducible to a single 3D equation. 

In Gross' spectator model |lj, Faddeev's type equations are deduced from the Bethe-Salpeter 
equation and the transition operator is written in terms of the two-body transition operators. The relative 
time variables are then eliminated by putting, in each three-body propagator, all the " offmassshallness" on 
the only fermion which interacts before and after. The Lorentz invariance-cluster separability requirement 
is satisfied by applying suitable Lorentz boosts on the two-body transition operators. 

The main difference between our approach and Gross' approach comes from the fact that we are 
(presently) interested by the (123) bound states only, to be computed principally by using a single 3D 
equation (such as the one from which Faddeev starts). Instead of working with the two-body transition 
operators, we work with the two-body potentials. We can however present our approach in terms of the 
two-body transition operators. Our approximation of each two-body transition operator is then unique 
and does not depend on the operators which come in front and behind it in the expansion of the Faddeev 
equations. This enables us to describe our model by a single 3D potential equation, by a set of three 
Faddeev equations, or by the expansion of the three-body transition operator in terms of the two-body 
ones. Another feature of our model is to not introduce Lorentz boosts by hand: the Lorentz-invariance / 
cluster separability requirement is exactly satisfied if we do not truncate our potentials (practically, this 
means that the satisfaction of this requirement and the approximation of the potentials can be improved 
together). Let us go back to the three- fermion Bethe-Salpeter equation and neglect for simplicity the 
three-body irreducible kernel K123 . Defining 

*12 = G°iG° 2 K 12 $, • • ■ $ = $ 12 + $ 23 + $31 , (140) 

we get 

(1 - G°G°/f 12 )$ = $ 23 + $31 (141) 
$ = (1 + G?G° T 12 ) ( $23 + $31 ) (142) 
$12 = G?G°T 12 ( $23 + $31). (143) 

Writing then 

$12 = G?G§G§ ft/fc/Wia, ••■ (144) 
in order to factor out the propagators, we get 

012 = ftft T u G?G° fcfa ( + 03i )■ (145) 

The (12) transition matrix element corresponding to our approximation (^4j) of Kyi is 

r i2(Pl20: P120, Po -P30) ~ T° 2 {p 30 ) 
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= K^+ + ( Sl2 , Sl2 ,P ^h 3 )(l-G 1 G 2 K^ ++ ( Sl2 , Sl2 ,P -h 3 ))- 1 . (146) 
This 2i 2 (p3o) is analytical in the Im(p3o) < half plane and 

T? 2 (h 3 ) = T+ + (s 12 ,s 12l P ~h 3 ). (147) 



This approximation, combined with eq uatio n ( [L45[) , implies that 0i2(pi2OjP3o) is independent of p\ 2 o 
(let us write thus (f>i 2 (p 3 o) )■ Equation (|l45| ) becomes then 

0i 2 (P3o) = /3i/%2? a (P3o) / dpi2oGSG°ft#! [ fetPio) + foi(pao)] (148) 



where pio,P20 must be written in terms of Po,P30:Pi20- Equation ( |148| ), together with similar equations 
for (f> 23 {pio) and 4>3i(P2o), admits solutions which are analytical in the lm(pko) < half planes. At 
Pm = hk, we get the 3D Faddeev equations 

<M/* 3 ) = T* 2 D++ (P - h a ) 1 . - [ ha(hi) + 3 i(M], ■ ■ ■ (149) 

r'o — o + ier 

which can easily be transformed back into our basic three-cluster potential equation (by performing in 
the reverse order the transformations made above at the 4D level). Note that we would get the same 
result by approaching directly T\ 2 by Ti 2 {h 3 ) instead of Tf 2 (f3o)- This means approaching K\ 2 by a 
given function K® 2 (j) 3 o), analytical in the Im(p3o) < half plane and equal to K^ 2 ++ (s\ 2 , s\ 2 , Po~h 3 ) 
at p 30 = h 3 . 



The key point of the manipulations above is the dominance of the positive-energy poles of G\G 



in (145). This was obtained by approaching K( 2 (or T\ 2 ) by a constant with t±[ 2 positive-energy 



projectors. We could try to insure the dominance of these two poles more economically. Let us come 



back to equation (145) without making any approximation on T\ 2 . The elements of (145) are then the 
operator and functions 

r i2 (K20! P120, Po -Pao) (150) 

012 (P'l20! P30 ), ^23 (P230, PlO ), 031 ( P310 , P2Q ) (151) 



and (145) must be integrated with respect to Pi 2 q, with p 3 Q fixed. We must thus write p 23 o, pio, paw, P20 
in terms of Pi 2 q, p 3 o- Searching as above for solutions <pij (p^ , Pko ) which are analytical in the lm(pko) < 
half planes, we shall close our integration path clockwise (counterclockwise) in front of 023 ( 03 1 ) and 
keep only the resi due of the pole of Af 2 + GlG 2 which puts fermion 1 (2) on its positive-energy mass shell. 



The elements of ( 145 ) are then replaced by 



T i2 (i>'i2o> P-lw, Po -P30) -> T 12 {p' 120 , S12 T P ° 2 ± P3 ° 2 h3 , Po - Pao) (152) 

^120 G°G° A/32 - ow A " + >,^- ( 153 ) 
(P - S) - {p 30 - h 3 ) +ie 

P — S 

ha (P230, Pio ) -> 023 ( S23 H 2 ( p30 ~ hz )' El ) ( 154 ) 

p> g 

hi (Paw, P20 ) -> 03i («3i ^ 1- (P30 -ha), E 2 ) (155) 



where we take the upper signs in (152) in front of (£> 23 , the lower signs in front of 4>3i- 

The manipulations above are submitted to some restrictions on the . In the p^o variable the Ty 
must be analytical in the lm(pko) < half planes. In the and pijo variables the Ty must be 

analytical in the whole complex plane. Moreover, the Ty must also be asymptotically bounded in the 
three variables and contain positive-energy projectors. These conditions are of course not satisfied 

by the exact transition matrix elements, but we shall assume that the singularities and the negative- 
energy parts of the can be neglected in the computation of the integrals with respect to the relative 
energies. 

Taking the equation for <f>i 2 at p 3 Q — E 3 and p' 12Q — s^ 2 — s\ 2 ± \{Po — S), means that we put 
fermions 2 and 3 (for +) or fermions 1 and 3 (for -) on their mass shells. This leads to 

4>\ 2 = -2iTT AA ( T^ 1 - 1 , 4>la + T? 2 2 1 . fa ) (156) 
Po — b + te Po — S + i€ 
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# 2 = -21TT (3^ ( T,Y 1 , celt + ^i2 2 p It ■ ) (157) 

where the upper indexes refer to the initial and final on mass shell non-spectator fcrmions. Equations 
( |l56| ),( [l57| ), together with similar equations for the (23) and (31) pairs, are a closed system of six 3D 
equations. In each Ty two fermions are on the mass shell: the spectator fermion k and the fermion 
which is not going to interact next left (p- j0 ) or next right (Pyo )■ This is of course the philosophy of 
Gross' spectator model [|l7j pi . 

If we replace s ± by s in the we get our basic three-cluster model. This supplementary approxi- 
mation can be justified by noticing that we already neglected the contributions of the singularities of Ty 
in the and pyo complex planes. The dependence on these variables can thus not be very strong, as 
we know that a function which is analytical and bounded on the whole complex sphere must necessarily 
be a constant. Using the same argument of the consistency of the approximations, we could also argue 
that Tij constant implies Kfj and constant with K^. = Kij. Our basic three-cluster equation is 
thus not a priori a worse approximation of the three-fermion Bethe-Salpeter equation than Gross' equa- 
tions, nor a better approximation than the simpler positive-energy instantaneous approximation or the 
corresponding Born approximation. It reflects the choice of preserving the exact cluster-separated limits. 

To conclude the comparison between our 3D reduction and Gross' spectator model, we can say that, 
at lowest order, our 3D reduction is much simpler: our two-body transition operators are always taken 
at the same value of the corresponding relative energie, so that we can transform the set of coupled 
equations for the transition operators into a unique equation for a wave function. We do not use Lorentz 
boosts, as we choose to exploit the implicit Lorentz covariance of the unapproximated separated-clusters 
limit equations, so that the effects of these boosts will be approached in principle by the progressive 
introduction of the higher-order two-fcrmion terms. Gross' higher-order corrections introduced by taking 
the transition operators with different fermions on the mass shell are not a priori more important than 
the contributions of the singularities of these transition operators themselves to the closed path integrals 
in the relative energies. Our brief calculation with the one-photon exchange graph (subsection 5.3) shows 
that, in this case, the contributions of these singularities are in fact more important, by an order of 
a (but it turns out that the leading order of these contributions will be cancelled by the leading order 
contribution of the first crossed graph (42) ). Our higher-order calculation takes both types of contributions 
into account. 



7 Heavy mass limits and external potentials. 

7.1 Heavy mass limits in the two- fermion problem. 

When the mass of one of the fermions goes to infinity, we must find the equation of the other fermion in 
a potential (a Dirac equation with a Coulomb potential in QED). This limit can be obtained directly, or 
via a rearrangement of the equation. With Breit's equation with the second-order constraint and with 
Sazdjian's equation, the higher-order ladder and crossed terms cancel mutually at the one-body limit, so 
that the correct limit of the potential is already contained in the Born term. With Salpeter's equation with 
the second-order constraint, and with the positive-energy equation with a second-order constraint, the 
limit is also a Dirac-Coulomb equation, but solved with respect to the A + ip part of the wave function: 
the A~ip part is transformed into higher-order contributions to the potential. The physical content 
remains identical to that of a Dirac-Coulomb equation, but the correct limit of the potential is no more 
entierely contained in the Born term, so that a truncation of the potential would spoil the one-body limit 
[8-12]. With the first-order constraint, the results are similar, but after complicated rearrangements of 
the equation. 

In the three-fermion problem, we expect similarly that, at the limit of an infinite fermion mass, we 
would get an equation describing two fermions in an external potential. It is thus interesting to compare 
this equation with the equation we wrote directly |]l6| for the two-fermion in an external potential system. 

7.2 Two fermions in an external potential. 

The two-fermion plus potential problem can be approached in two ways which we shall call the two- 
body and the three-body approaches. In the two-body approach, the external potential is included in the 
definition of the creation and annihilation operators of field theory. Practically, we can keep the equations 
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obtained in the treatment of the pure two-fermion problem, using simply a generalized definition of the 
hi 0,0,0: 

hi = di .'pi + A Tii + Vi(xi) (158) 

where is the external potential acting on the fcrmion i. All quantities can be expanded on the basis 
built with the eigenstates of hi and hi- In the three-body approach, the free creation and annihilation 
operators are used and the external potential is treated as an heavy third body. This approach will be 
presented as an heavy mass limit of the three fermion problem in the next subsection. 

The equations for two fermions in an external potential must be written in the laboratory reference 
frame (in which the external potential is defined) and are not Lorentz invariant. If we "switch off' the 
mutual interaction we get a pair of uncoupled equations for two independent fermions in an external 
potential. If we switch off the external potential the equations remain written in the laboratory frame, 
which still refers to the vanished external potential. If the reduction series is not truncated, this no- 
external interaction limit of the equation is equivalent to a pure two-fermion covariant Bethe-Salpeter 
equation. 



7.3 Heavy mass limits in the three-fermion problem. 



Let us now take our 3D equation (p2[), with the potentials Vij obtained by performing 3D reductions 
based on the G^ p as in section 2, but with the second-order constraints <5(pyo — Mij) instead of the 
first-order constraints <5(pyo — Sij)- When 7713, for example, goes to infinity (two-body limit), we get, 
writing P — W\% + m 3 : 

* = w — \ — r ( yi2 + y 2 + + i> (159) 

w 12 — Hi — n% 

with (Po — /13) replaced by Wi 2 in V\ 2 - The potential V 2 + is given by the series 



A 



V 2 + V 2 



{W 12 -h 1 ) + E 2 



V 2 



A 2 + 



(160) 



where V% is an external potential acting on fermion 2 (it is a Coulomb potential in QED, if we use the 
second-order constraint, but in general it could still depend on (W12 — hi)). The potential V^~ is given 
by a similar formula. 

The differences with the two-body approach of the two-body plus potential problem are: 

- The projectors A^ are the free ones. 

- The external potential for fermion 2 is now V 2 + , given by the series (16C). This is the potential 
obtained by solving the equation for fermion 2 in the external potential V 2 with respect to the positive 
free energy components, the unobservable energy of fermion 2 being replaced in it by W\ 2 — h\. 

-V 2 itself could still depend on (Wi 2 — h\) and vice-versa. 



The cluster separability property in the three-body way survives the high-mass limit: switching off V 2 
and Vi leads to the equation for the (12) two-fermion system, switching off Vi 2 and V\ leads to a free 
Dirac equation for fermion 1 with the A 2 projection of the Dirac equation for fermion 2 in the external 
potential V 2 . In the two-body approach of the two-fermion plus potential problem we had however more: 
switching off the mutual interaction only led to a pair of independent equations for each fermion in the 
external potential. Here, in the three-body approach, the equation does not split perfectly into two 
parts: we can write W\ 2 = W\ + W 2 , but we have W\ 2 — hi instead of W 2 in the serie s defining V 2 
and vice- versa. This is a consequence of the energy dependence introduced into (160) by the use of 
the anti-continuum dissolution projectors, combined with our choice of putting the spectator fermion on 
the mass shell in each two-body interaction, and neglecting the three-body terms which could balance 
this modification. The discrepancy is of order V 4 . At the two-fermion plus potential level, a better but 
still not perfect equation would be obtained by replacing hi by hi + Vi in the series defining V 2 and 
vice- versa. 



8 Conclusions 

We have written a 3D equation for three fermions (our basic three-cluster equation) by combining the three 
two-body potentials obtained by an untruncated propagator-based 3D reduction of the corresponding 
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two-fermion Bethe-Salpeter equations, putting the spectator fermion on the mass shell. In this way, 
the cluster-separated limits are exact, and the Lorentz invariance / cluster separability requirement is 
automatically satisfied, provided no supplementary approximation, like the Born approximation, is made. 
Such a truncation would render the Lorentz covariance of the two-fermion clusters only approximate (a 
Born approximation preserving the Lorentz invariance / cluster separability property can be obtained by 
using another 3D reduction based on a covariant second-order two-body propagator of Sazdjian, combined 
with a covariant substitute of A ++ . This leads to a 3D three-cluster equation which is covariantly Born 
approximable, but more complicated |4p| ). Our 3D equation can be written in terms of the two-body 
potentials, or in terms of two-body transition operators (Faddcev formalism). The use of positive free- 
energy projectors in the chosen reductions of the two-fermion Bethe-Salpeter equations prevents our 
3D three-fermion equation from continuum dissolution. The potentials are hermitian and depend only 
slowly on the total three-fermion energy. The one high-mass limits of our "basic three-cluster equation" 
are approximately exact. The correction of the remaining discrepancy would demand the introduction of 
higher-order three-body terms. 

Our combination of cluster separability and Lorentz invariance in the three fermion problem makes 
explicit use of the fact that the clusters can only be two-fermion and/or free fermion states. This is not 
directly adaptable to four or more fermion systems. In this respect, 3 is still not N. 

Our 3D equation can be written directly, or derived from an approximation of the three-fermion 
Bethe-Salpeter equation, in which the three-body kernel is neglected while the two-body kernels are 
approached by positive-energy instantaneous expressions, with the spectator fermion on the mass shell. 
The correction terms are thus known at the Bethe-Salpeter level and can be transformed into corrections 
to the 3D equation. 

Our lowest-order three-fermion 3D equation is much simpler than the equations written in the frame- 
work of Gross' spectator model, which contain some induced higher-order three-fermion contributions. 
Our calculation of the next order shows that these contributions are not a priori more important that the 
ones coming from the singularities of the kernels, neglected in Gross' model. We think that the relative 
importance of the different contributions could largely depend on the specific problem studied. 

There exists an infinity of ways of performing the 3D reduction of the two-fermion Bethe-Salpeter 
equation. The potentials generated by these reductions could all by used to build a three-fermion 3D 
equation, keeping however in mind the continuum dissolution problem. This leaves us a large freedom to 
suit phenomenological needs. 

We have also examined the possibility of using kernel-based 3D reductions in the two and three- 
fermion problems. We found that this kind of 3D reduction can also be used in the calculation of bound 
states, despite the fact that they do not lead to the correct scattering operators (excepting in Phillips 
and Wallace's reduction). 

As a first test of our methods in a specific problem, we are presently working on a system of two or 
three fermions exchanging photons in Feynman or Coulomb's gauges [fi"2"|| . The various terms generated by 
our propagator-based 3D reduction will be compared with these generated by a kernel-based 3D reduction 
and with these written in the framework of Gross' spectator model. 
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